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基于 离散 型 伴随 方法 的 单 边 膨胀 喷 管 优化 设计 研 
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宋 红 超 FÈ 季 路 成 

(北京 理工 大 学 ， 北 京 ，100081) 
摘要 : 目前 高 超声 速 喷 管 的 设计 方法 主要 有 特征 线 方法 与 基于 CFD 分 析 的 优化 设计 方法 。 采 用 特征 线 法 设计 喷 管 时 流动 需 
基于 无 粘 、 无 旋 假 设 进行 简化 ,无 法 确保 设计 结果 在 真实 流动 情况 下 达到 最 优 : 采用 遗传 算法 或 常规 梯度 类 方法 等 优化 算 
法 进行 优化 设计 时 ， 设 计 变 量 数目 的 增多 为 设计 带 来 计算 量 方面 极 大 的 挑战 。 为 克服 上 述 方法 在 设计 过 程 中 的 不 足 ， 依 托 
伴随 方法 具有 计算 量 与 设计 变量 数目 几乎 无 关 的 特点 ,本 文 发 展 了 基于 离散 型 伴随 方法 的 单 边 膨胀 喷 管 的 设计 方法 ， 以 控 
制 喷 管 沿 流程 的 面积 分 布 作为 参数 化 方法 , 在 初步 设计 的 性 能 优良 的 原型 喷 管 基础 上 实现 了 精细 化 伴随 优化 设计 , 优化 后 
推力 系数 比 原型 喷 管 提升 0.8 个 百分点 。 
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Research on the Optimization of Unilateral Expansion Nozzle Based on the 


Discrete Adjoint Method 


HongChao Song! Li Xin? Ji Lucheng? 
(Beijing Institute of Technology, Beijing, 10081, China) 

Abstract: At present, the characteristic line method and the optimization method based on the CFD analysis are 
the mainly methods for hypersonic nozzles designing. However, the characteristic line method can’t come over 
the shortcoming of the assumptions that the flow is no-viscosity and no-rotation. Besides, the genetic algorithm 
and general gradient-based algorithms used for hypersonic nozzles designing can’t handle the challenges of 
increasing number of design variables. Considering the calculation cost of the adjoint method isn’t sensitive 
with the number of design variables, this paper uses the discrete adjoint method to design the unilateral 
expansion nozzle, accomplishing the parametric model of controlling the area distribution of the nozzle along 
the streamline. On the base of the initial nozzle designed by the genetic algorithm, which already has the great 
performance, the optimization can improve the thrust coefficient with 0.8 percentage points. 

Key Words: Discrete Adjoint Method; Unilateral Expansion Nozzle; Parametric Model 


0 引 设计 方法 , 该 方法 在 给 定 喷 管 长 度 和 质量 流量 的 约 
束 条 件 下 ,采用 特征 线 法 计算 流 场 ,， 并 据 此 获得 最 


单 边 膨胀 喷 管 广泛 应 用 于 高 超声 速 飞行 器 种。 大 推力 型 面 。Zudov 等 构造 并 求解 了 以 超声 速 一 
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压 发 动机 的 喷 管 中 。 目前, 针对 此 类 高 超声 速 喷 管 维 非 对 称 喷 管 最 大 推力 为 目标 函数 的 变 分 问题 ,其 


mi YA MEI de^ 4x RE ER 而 AKI 
的 设计 方法 主要 包括 特征 线 法 以 及 优化 设计 方法 。 — 求解 过 程 使 用 了 特征 线 万 法 。 国 内 ， 曹 德 一 等 
也 采用 特征 线 方法 直接 对 后 体 单 边 膨胀 喷 管 进行 
特征 线 法 被 广泛 用 于 超声 速 流动 型 面 设计 中 ， 。 了 设计 。 
可 设计 出 内 部 流 场 光滑 且 不 存在 激 波 的 喷 管 .特征 | 
线 法 的 一 种 典型 方法 是 Rao"" 提 出 的 最 大 推力 喷 管 早期 特征 线 法 在 高 超声 速 内 流 部 件 型 面 设计 
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中 取得 一 定 效果 。 然 而 , 特征 线 法 需要 基于 流动 无 
粘 、 无 旋 的 假设 , 与 真实 流动 具有 较 大 的 偏差 , 同 
时 利用 特征 线 法 设计 出 的 喷 管 通常 情况 下 管道 较 
长 ， 需 要 对 设计 结果 进行 截 短 ,导致 设计 喷 管 具有 
较 大 的 性 能 损失 。 


随 着 计算 流体 力学 和 计算 机 技术 的 迅猛 发 展 ， 
近 几 十 年 出 现 的 基于 CFD 分 析 的 优化 设计 方法 取 
得 了 重大 进步 。 陈 兵 “ 基 于 NS 方程 的 优化 设计 建 
立 了 二 维 尾 喷 管 模型 ,应 用 复合 形 法 对 该 模型 进行 
了 求解 , 随后 又 结合 遗传 算法 与 高 效 、 高 精度 的 空 
进 算法 “， 对 二 维 超声 速 超 燃 冲压 发 动机 尾 
行 了 气动 优化 设计 研究 。 甘 文 彪 等 ” 则 结合 
设计 方法 、 痊 代 模 型 技术 以 及 遗传 算法 构建 了 
优化 方法 , 并 将 其 应 用 于 高 超声 速 飞行 器 后 体 
喷 管 的 一 体 化 设计 中 。 
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上 述 高 超声 速 喷 管 设计 中 所 采用 的 遗传 算法 
以 及 常规 梯度 类 方法 等 , 尽管 有 效 , 但 计算 量 均 会 
iioi inca de DA EA pd 增长 ， 
难以 满足 全 三 维 工 程 设计 中 缩短 设计 周期 、 提 高 设 
计 自 动 化 以 及 精细 化 的 需求 。 


相 比 常规 优化 方法 , 伴随 方法 具有 计算 量 与 设 
计 变 量 数目 几乎 无 关 的 明显 优势 , 能 够 克服 常规 优 
化 方法 在 计算 量 方面 的 不 足 。 在 Jameson 和 Giles 
等 人 “的 努力 下 ， 经 过 二 十 多 年 的 发 展 ， 辟 型 、 
机 翼 、 嗓 身 组 合体 乃至 全 机 的 气动 伴随 优化 设计 已 
得 到 较为 成 熟 的 应 用 。 本 文采 用 离散 型 伴随 优化 方 
法 设计 高 超声 速 单 边 膨胀 喷 管 , 真正 以 精细 化 设计 
模式 全 面 解决 高 超声 速 单 边 膨胀 喷 管 设 计 问题 ,并 
为 后 续 的 进 气 道 和 推进 系统 一 体 化 设计 过 程 提 供 
了 技术 参考 。 
1 离散 型 伴随 优化 方法 
1.1 离散 型 伴随 优化 理论 

对 于 定常 流动 优化 问题 ， 流 动 方 程 的 残 差 R 
可 视 为 优化 设计 变量 D 、 网 格 X 和 流动 变量 Q 的 
函数 ， 即 : R- R(DQ,X)- 0。 上 述 的 残 差 方程 
本 质 上 属于 定常 流动 优化 问题 的 限制 条 件 。 在 初始 
目标 函数 A(D.Q, X) 中 引入 拉 格 朗 日 乘 子 / ， 新 目 
标 函 数 工 为 : 


As. d /^EHRTI| 
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L(D,Q,X,4)- f(D,Q,X)+ XR(D,Q,X) 
(1) 
x (1) 将 约束 优化 问题 转化 为 无 约束 优化 问 
题 ， 其 中 f(D, O, X) 为 气动 目标 函数 (一 般 求 最 小 


值 ) 1 是 拉 格 朗 日 乘 子 矢量 (伴随 变量 ) 。 式 (1) 
对 DD 求 导 : 


pe a 
pii as Px? 


Euh: Jb NU RU HER 


W 


大 的 主要 原因 。 由 于 拉 格 朗 日 乘 子 的 任意 性 ， 
REMOS 的 系数 为 0 可 以 相应 地 简化 计算 


二 入 = - — 3 
bo! Q = 


方程 (3) 即 为 流动 问题 的 离散 型 伴随 方程 。 满 
足 伴随 方程 的 基础 上 ， 简 化 后 的 目标 函数 为 : 


IN dos 

nox 

B (4) 
PE or pu 

1.2 流动 约束 方程 


1.2.1 流动 控制 方程 组 
维 无 量 纲 雷 诺 平 均 NS 方程 可 以 表述 为 如 下 


(K 


Bp: 


oQ 
Vo +f a mar -if c mar = 0 (5) 


其 中 , n 为 控制 体 边 界 的 外 法 向 单位 向 量 ，Q 
为 单元 平均 守恒 变量 ， 下 为 无 粘 通 量 ， 耻 , 为 粘性 
通 量 。 层 流 粘 性 计算 采用 完全 气体 假设 , 求解 遵循 
Sutherland 公式 。 
1.2. 2 满 流 模型 


经 由 大 量 实验 表明 ， 本 文选 用 的 SA ini E78 
具有 非常 强 的 鲁 棒 性 , 极 适 用 于 一 般 内 流 航空 工程 


的 应 用 。 无 量 纲 SA 滑 流 模型 控制 方程 如 下 : 
p2 = L (V [0 (16, )0V 9] 
~ M C, D 
—c, OV? D — —&. 一 一 人 6 
ao- f -m (T) — € 
+c, (1- f, )S0 
其 中 : 
Ü 
d 之 ?N= 一 
PT D 
" M, 
NE E 二 > Ja cde E 
14 cf 
X 5 
=1 
^ l4 xf, f T 十 C | 
M 加 
二 了 十 S s — LEN MN 
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(7) 


式 (6) 中 5 为 涡 量 的 模 ，4 为 当前 网 格 单元 到 
壁面 的 最 短 距离 。 
1.3 优化 方法 

本 次 工作 开发 的 优化 程序 包括 问题 目标 函数 建立 、 
网 格 生 成 、 流 场 求解 、 伴 随 场 求解 、 敏 感性 计算 、 优 化 
算法 、 网 格 变形 共 七 个 模块 ， 详 细 优 化 流程 (如 流 场 求解 
方法 、 伴 随 场 求解 方法 、 优 化 算法 以 及 网 格 变形 等 ) 参见 
文献 12。 
1.3.1 目标 函数 


为 使 高 超声 速 喷 管 推力 系数 最 大 化 ,本 文选 取 
的 目标 函数 形式 为 : 


f--TIT, (8) 
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管 为 超声 速 进口 , 依 特征 线 理论 , 进口 需 指 
定 所 有 的 气动 参数 。 进 口 给 定 马 赫 数 1.4， 静 压 
145600 Pa ， 静 温 1500K ， 比 热 比 1.4; 出 口 选择 
压力 边界 , 并 假设 出 口 边界 绝热 和 等 粹 , 如果 出 口 
为 超声 速 流动 , 则 所 有 的 气动 参数 均 采用 外 推 的 方 
式 进行 计算 ; 壁面 条 件 采 用 无 滑 移 壁面 。 当 前 程 请 
初始 条 件 选 择 以 来 流 参数 对 全 场 进 行 初始 化 赋值 。 
1.3.3 几何 参数 化 方法 

本 文 的 设计 变量 选取 喷 管 表面 的 500 个 节点 ， 
当 发 生 扰动 时 , 每 个 设计 变量 都 在 垂直 于 流 线 方向 
的 平面 内 发 生变 动 , 继而 改变 喷 管 沿 流程 的 面积 2 
布 。 因此 , 参数 化 对 象 为 气动 形状 的 扰动 而 非 气动 
形状 本 身 , 通过 控制 喷 管 的 H 
Ia aae. T NER. nh 
管 参数 化 方式 为 : 


nV)= r+ Dr,Q) (9) 

其 中 , vy 为 包含 设计 变量 的 矢量 ，x,(v) 是 喷 管 

型 面 在 点 n ABA ES. r^ 为 初始 模型 的 型 面 坐标 ， 
Dr,(v) 是 扰动 坐标 ， 扰 动 坐标 的 具体 表达 式 为 : 


Dr()- à vn? (10) 
il 


其 中 , m 表示 设计 变量 的 个 数 ，p,, 表示 第 i 个 
设计 变量 的 指数 ，n, 表示 在 点 n 处 的 单位 法 向 矢 


=E. 
FH o 


2 原型 喷 管 方案 设计 


梯度 类 算法 属于 局 域 优 化 方法 的 一 种 , 其 最 终 
优化 目标 为 确定 目标 函数 在 优化 域内 的 一 个 极 值 
而 非 最 值 , 不 能 保证 优化 结果 的 全 局 最 优 性 。 针 对 
此 问题 ,一 种 可 行 的 解决 方式 为 预先 确定 最 优 解 所 
在 的 局 部 优化 域 , 缩小 伴随 优化 的 优化 区 间 以 保证 
优化 结果 的 有 效 性 ， 因 此 原型 喷 管 的 设计 至 关 重 


推力 系数 , 也 是 初始 喷 管 的 推力 


系数 。 
1.3.2 初 边界 条 件 


本 次 设计 的 单 边 膨胀 喷 管 以 4. 5 马赫 的 速度 


要 。 原 型 喷 管 作为 伴随 优化 的 起 点 , 起 到 了 “基因 
RERE” REH, 其 性 能 一 定 程度 上 决定 了 后 续 伴 
随 优化 取得 的 效果 。 


飞行 在 18. 5km 的 高 空 ， 查 得 此 高 度 的 标准 大 气 吏 
压 为 6994. 8Pa。 


为 了 提高 性 能 , 原型 喷 管 的 设计 应 采用 具有 全 
RSEN a 为 此 , 控制 原始 喷 管 造 型 


要 尽 可 能 少 以 减 小 原型 喷 管 的 设计 计 香 
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量 。 本 文 的 原型 喷 管 设计 过 程 中 , 给 定 尾 喷 管 进口 
截面 为 圆 形 , 喷 管 出 口 为 矩形 , 整个 优化 过 程 中 限 
制 进出 口 形状 以 及 喷 管 长 度 不 变 。 在 此 基础 上 , 沿 
流向 的 各 展 向 剖面 取 为 超 椭圆 , 其 面积 和 mn 指数 沿 
流程 分 布 作为 优化 参数 , 优化 过 程 中 , 通过 调节 面 
积分 布 规律 、 超 椭圆 指数 n 分 布 规律 以 及 下 底板 角 
度 等 9 个 参数 ， 即 可 唯一 确定 喷 管 形状 ， 见 图 1。 


图 1 喷 管 造型 示意 图 


Fig. 1 Diagram of Nozzle Modeling 


VETE PIUSUORE gg D E bevr P I — HER TUS 
制 能 力 ,特别 注重 面积 沿 流向 的 分 布 规律 以 及 二 面 
角 原 理 的 应 用 。 

采用 遗传 算法 进行 原型 喷 管 的 设计 时 , EDS RI] 
网 格 数量 为 60 万 左右 ， 为 保证 计算 精度 ，Y-plus 
的 值 小 于 1， 网 格 见 图 2。 


(a) 壁面 网 格 布置 (b) 进出 口 网 格 布置 


图 2 喷 管 网 格 布置 


Fig.2 Nozzle Grids 


在 9 个 控制 参数 的 基础 上 , 本 文 利用 遗传 算法 
进行 进化 寻 优 使 喷 管 原始 造型 接近 全 局 最 优 解 . 遗 
传 算法 的 种 群 规 模 为 10。 由 图 3 可 发 现 种 群 进化 
到 第 八代 时 ， 推 力 系数 由 不 足 91% 迅 速 提升 至 93% 
左右 , 且 优 化 过 程 趋 于 稳定 。 将 遗传 算法 优化 过 程 
的 第 八代 计算 结果 作为 原型 喷 管 的 最 终 设 计 , 为 后 
续 伴 随 优化 提供 原始 造型 .图 4 为 原型 喷 管 马赫 数 
在 流动 对 称 面 内 的 分 布 云图 以 及 治 流程 的 分 布 云 
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图 ， 可 以 发 现 ， 出 口 处 马赫 数 分 布 较为 均匀 。 


9904 4 5 
优化 次 数 


图 3 遗传 算法 收敛 史 


Fig.3 Convergence History of the Genetic 


Algorithm 


(b) 马赫 数 沿 流程 的 分 布 


图 4 原型 喷 管 流 场 马赫 数 云图 


Fig. 4 Mach Number Distribution of the Original Nozzle 
3 伴随 优化 结果 分 析 


以 原型 喷 管 为 起 点 , 通过 喷 管 表面 的 500 个 设 
计 变 量 对 喷 管 进行 精细 化 伴随 优化 设计 ,这 种 精细 
化 程度 是 遗传 算法 及 其 它 常 规 梯度 类 算法 难以 实 
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现 的 。 图 5 为 目标 函数 随 优化 过 程 的 迭代 曲线 。 可 
以 发 现 , 在 前 三 个 优化 迭代 过 程 中 , 目标 函数 的 提 
升 较 慢 , 这 与 优化 算法 、 目 标 函 数 形式 以 及 参数 化 
方法 都 有 关系 ; 同时 , 优化 在 第 四 个 优化 过 程 中 取 
得 了 最 大 的 效果 ， 推 力 提升 了 0. 8 个 百分点 左右 ; 
而 后 目标 函数 基本 收敛 , 这 预示 着 优化 已 经 达到 了 
极 值 ， 优 化 过 程 结束 。 本 次 伴随 优化 只 进行 了 7 
次 优化 循环 ,计算 量 仅 相当 于 14 倍 流 场 的 计算 量 。 


4 
优化 次 数 


图 5 伴随 优化 目标 函数 收敛 史 


Fig.5 Convergence History of the Objective Function 


Optimized by the Adjoint Method 


pd 


优化 前 


(a) 对 称 面 对 比 
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(b) 二 面 角 区 对 比 


6 优化 前 后 几何 模型 对 比 


Fig.6 Comparison of Geometry before and after 


Optimization 


本 次 伴随 优化 过 程 采 用 了 常规 的 几何 约束 方 
法 ， 限 制 进出 口 几何 保持 不 变 。 图 6 为 优化 前 后 
几何 模型 对 比 , 可 以 发 现 , 优化 前 后 几何 模型 变化 
非常 大 ， 其 中 ,圆圈 处 标记 了 变化 最 大 位 置 处 。 图 
6 (b) 显示 了 二 面 角 区 的 变化 尤为 明显 。 


(a) 原型 喷 管 对 称 面 马赫 数 云图 
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(b) 优化 后 对 称 面 马赫 数 云图 


(c) 原型 喷 管 对 称 面 静 压 云图 


(d) 优化 后 对 称 面 静 压 去 


R] 


7 优化 前 后 对 称 面 云图 对 比 


Fig.7 Comparison of Flow Field before and after 


Optimization 


图 7 为 优化 前 后 流动 对 称 面 内 马赫 数 和 静 压 
分 布 云图 。 可 以 发 现 , 优化 前 后 的 流 场 变化 比较 明 
显 : 原型 喷 管 从 进口 到 出 口 几 何 变化 平缓 , 流 场 也 
是 比较 平缓 地 膨胀 ; 相 较 于 优化 前 , 优化 后 喷 管 前 
半 部 分 流 场 膨 胀 程度 较 低 ， 而 后 半 部 分 则 急剧 膨 
HE. 下 底板 表现 尤为 明显 。 图 8 展示 了 马赫 数 沿 流 
程 的 分 布 , 也 可 以 反映 出 优化 后 喷 管 前 半 部 分 膨胀 
程度 小 , 但 到 出 口 处 , 优化 后 喷 管 的 膨胀 程度 反而 
比 原 型 喷 管 更 大 ; 尤其 在 角 区 附近 , 优化 后 的 喷 管 
出 口 四 个 角 区 膨胀 程度 更 充分 , 角 区 流动 得 到 了 极 
大 改善 。 
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b) 优化 后 马赫 数 沿 流程 的 分 布 


对 


8 优化 前 后 马赫 数 沿 流程 的 分 布 


Fig.8 Comparison of Mach Number Distribution along 


the Streamline before and after Optimization 


4 总 结 


本 文 首先 用 较 少 的 参数 构造 喷 管 ,利用 遗传 算 
法 进行 全 局 性 优化 设计 ， 得 到 性 能 良好 的 原型 喷 
管 ; 继而 采用 离散 型 伴随 方法 对 原型 喷 管 进行 精 旨 
化 优化 设计 ,得 到 了 推力 系数 继续 提升 0.8 个 百 分 
点 的 最 终 设计 结果 。 


本 文 开 发 的 利用 伴随 方法 对 单 边 膨胀 喷 管 进 
行 设 计 的 体系 , 无 需 对 流动 进行 简化 假设 , 并 且 不 
受到 设计 变量 数量 的 限制 , 克服 了 现 有 设计 方法 的 
部 分 缺陷 ， 满 足 了 精细 化 设计 的 发 展 趋势 。 同 时 ， 
本 文 的 参数 化 模型 为 进 气 道 及 推进 系统 一 体 化 优 
化 设计 提供 了 关键 技术 思路 。 


参考 文献 : 


[1]Rao G V. Exhaust Nozzle Contour for Optimum Thrust[J]. 


Jet Propulsion: 1958,V28:377-382 

[2]Zudov V N,Lokotko A V ,Rylov A LNumerical and 
Experimental Investigation of Two-Dimentional Asymmetric 
Nozzles[R].AIAA: 96-3141,1996 

[3] 曹 德 一 ， 李 椿 营 . 高 超声 速 飞行 器 尾 喷 管 的 优化 设计 加 . 
北京 航空 航天 大 学 学 报 : 2007，33(10): 1162-1165 

Cao Deyi Li ChunXuan. The Optimation and Design of 


Outlet Nozzle in Hypersonic Aerocraft[J]. Journal of Beijing 
University of Aeronautics and Astronotics: 2007, 33(10): 
1162-1163 

[4] Fc, IRIE, 3 EL. — PRGER A LH e s LEE E Ai 
[J]JESETCR:2002, 23(5) 433-437 

Chen Bing, Xu Xu, Cai Guobiao. Optimization And Design of 


Tail Nozzle for two dimensional scramjet[J].Journal of 
Propulsion Technology: 2002,23(5): 443-437 

[5] BRES, B, 3 ERU. ek FAE SETA t E E338 77 TS BI] E B 
扩张 喷 管 优 化 设计 研究 中 .航空 学 报 :2007，28(4):827-832 
Chen Bing, Xu Xu, Cai Guobiao. Optimization Design of 


Unilateral Expansion Nozzle Based on Genetic Algorithm and 
Space Propulsion Method[J]. Acta Aeronautica Astronautica 
Sinica: 2007,28(4):827-832 

[6] E 3C, EUER E ERE ET Je V Fé E UL VEU. 
北京 航空 航天 大 学 学 报 :2011，37(6): 1-6 

Gan Wenbiao, Yan Chao. Optimum Design of Hose / Tail 


Nozzle for Hypersonic Vehicle[J]. Journal of Beijing 
University of Aeronautics and Astronotics:2011, 37(6):1-6 
[7]Jameson A.Aerodynamic Shape Optimization Using the 
Adjoint Method[R]. Lectures: Brussels, 2003 

[8]Giles M B, Pierce N A. Adjoint Equations in CFD: Duality, 
Boundary Condition and Solution Behavior[J]. AIAA Paper: 
977-1850, 1997 

[9]Reuther J, Jameson A. Control Based Airfoil Design using 
the Euler Equations[J]. AIAA Paper: 94-4272-CP,1994 
[10]Reuther J, Jameson A.Aerodynamic Shape Optimization 
of Wing and Wing-body Configurations using Control 
Theory[J]. AIAA Paper: 95-0123-CP,1995 

[11]Reuther J, et al. Constrained Multipoint Aerodynamic 
Shape Optimization using an Adjoint Formulation and Parallel 
Computers[J]. AIAA Paper: 97-0103 

[12] 宋 红 超 ， 季 路 成 .目标 函数 形式 对 伴随 优化 影响 的 研究 . 
工程 热 物理 会 议 :2016. 


Song Hongchao, Ji Lucheng. Influences of the Form of 


Objective Function on the Result Optimizied by Adjoint 
Method[J]. Meeting Paper: 2016 


